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In this paper we examine ^-derivable approximations in QED. General theorems tell us that the gauge 
dependence of the n-loop <l?-derivable approximation shows up at order g^" where g is the coupling constant. 
We consider the gauge dependence of the two-loop ^-derivable approximation to the Debye mass and show 
that it is of order e'^ as expected. We solve the three-loop ^-derivable approximation in QED by expanding 
sum-integrals in powers of e^ and m/T, where m is the Debye mass which satisfies a variational gap equation. 
The results for the pressure and the Debye mass are accurate to order e^. 
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I. INTRODUCTION 

The thermodynamic functions for hot field theories 
can be calculated as a power series in the coupling con- 
stant g at weak coupling. The free energy has been cal- 
culated through to order g'^ in 1, 2] for scalar 0'* theory, 
in [3 for QED and in ^J for nonabelian gauge theories. 
The corresponding calculations to order g^ were carried 
out in Refs. 4,^, Refs. 0,01 and Refs. [alalj respec- 
tively. In Fig.^ we show the successive perturbative ap- 
proximations to P/7^idcai as a function of e(27rT). Each 
partial sum is shown as a band obtained by varying the 
renormalization scale /i by a factor of two around the 
central value /i = 2ttT. We have only done this varia- 
tion for the e^ approximation. To express e(/Lt) in terms 
of e(27rT), we use the solution to the one-loop renormal- 
ization group equation in QED. The Figure shows that 
the weak-coupling expansion is poorly convergent unless 
the coupling constant is small and that it is very sensi- 
tive to the renormalization scale /i. The lack of conver- 
gence seems to be related with screening and quasipar- 
ticles which is associated with the soft momentum scale 
of order eT. The instability of the weak-coupling expan- 
sion is a generic problem in hot field theories and makes 
it essentially useless for quantitative predictions. 

There are several ways of systematically reorganizing 
the perturbative expansion to improve its convergence 
properties and various approaches have been discussed in 
detail in the review papers Refs. [lOllllLll^ . One of these 
methods is screened perturbation theory (SPT) which in 
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FIG. 1: Weak-coupling expansion for the pressure to orders 
e^ (dotted curve), e^ (dashed curve), e* (long-dashed curve), 
and e^ (solid lines-|-band) normalized to that of an ideal gas 
as a function of e{2nT). 



the context of thermal field theory was introduced by 
Karsch, Patkos and Petreczky. 13] (See also refs. [liLllSl 
llQj). In this approach, one introduces a single variational 
parameter which has a simple interpretation as a thermal 
mass. In SPT a mass term is added to and subtracted 
from the scalar Lagrangian with the added piece kept 
as part of the free Lagrangian and the subtracted piece 
associated with the interactions. The mass parameter 
satisfies a variational equation which is obtained by the 
principle of minimal sensitivity. 

In gauge theories, one cannot simply add and sub- 
tract a local mass term as this would violate gauge in- 
variance. Instead one adds and subtracts to the La- 
grangian a hard thermal loop (HTL) improvement term. 



The free part of the Lagrangian then includes the HTL 
self-energies and the remaining terms are treated as per- 
turbations. Hard thermal loop perturbation theory is a 
manifestly gauge invariant approach that can be applied 
to static as well as dynamic quantities. SPT and HTL 
perturbation theory have been applied to three and two 
loops [Hill El IS mill, respectively, and the con- 
vergence properties are improved dramatically compared 
to the weak-coupling expansion. 

The ^-derivable approach is another way of reorga- 
nizing the perturbative expansion which is variational of 
nature. In this approach, one uses the exact propagator 
as a variational function. Its formulation was first con- 
structed by Luttinger and Ward ,2^] and by Baym [2J|. 
Later it was generalized to relativistic field theories by 
Cornwall, Jackiw and Tomboulis |2a|. The approach is 
based on the fact that the thermodynamic potential can 
be expressed in terms of the two-particle irreducible (2PI) 
effective action which has a diagrammatic expansion in- 
volving the 2PI skeleton graphs. Although here we focus 
on equilbrium physics we note that the 2PI formalism 
and its generalizations are also very useful when study- 
ing non-equilibrium real-time physics |2fil l27l l28l | . 

The ^-derivable approach has several attractive fea- 
tures. One is that it respects the global symmetries of 
the theory. Thus it is consistent with the conservation 
laws that follow from the Noether's theorem. Second, 
when evaluated at the stationary point, one is guaran- 
teed thermodynamic consistency |24| . Finally, it turns 
out that the two-loop ^-derivable approximation has an 
additional property. The entropy reduces to the one- 
loop expression at the variational point. This property 
was first shown for QED by Vanderheyden and Baym [23 
and later generalized to QCD by Blaizot, lancu and Reb- 
han I2IIM El- 
Applying the ^-derivable approach to quantum field 
theories, one is facing two nontrivial issues. The first is- 
sue is the question of renormalization. The three-loop 
calculations by Braaten and Petitgirard ^] for a mass- 
less scalar field theory indicate that there are ultraviolet 
divergences at order g^ that cannot be eliminated by any 
renormalization of the coupling constant. These calcu- 
lations seem to contradict the results from the papers 
by van Hees and Knoll |33, and by Blaizot, lancu and 
Reinosa |3a|, which show that the 2PI effective action 
can be systematically renormalized. 

The second issue is that of gauge-fixing dependence. 
While the exact 2PI effective action is gauge independent 
at the stationary point, this property is often lost in ap- 
proximations. The problem has recently been examined 
by Arrizabalaga and Smit |33, who showed that the n- 
loop (f>-derivable approximation $„, which is defined by 
the truncation of the action functional after n loops, has 
a gauge dependence that shows up at order 5^". Fur- 
thermore, if the nth order solution to the gap equation is 
used to evaluate the complete effective action, the gauge 
dependence first shows up at order g"*". For a general 



proof of the gauge invariance of the exact 2PI effective 
action we refer the reader to Ref. (3^] . 

In gauge theories, approximate solutions to the 
gap equations in two-loop <I>-derivable approximation in 
terms of HTL self-e nerg ies have been obtained by Blaizot, 
lancu and Rebhan |33,|3llll2l) and by Peshier [33. For 
hard external momentum, these solutions are obtained 
by evaluating the one-loop self-energy diagrams with bare 
propagators. For soft external momentum, the solution is 
simply the HTL self-energies, which is in the imaginary- 
time formalism reduces to the Debye mass. The resulting 
approximation reproduces the pressure to order g'^ . 

Finally, we mention that a dimensionally reduced 
version of the three-loop ^-derivable approximation was 
solved in the case of scalar field theory [l3|. This was 
done by treating the contribution to the pressure from 
the nonzero Matsubara modes in strict perturbation the- 
ory and applying the <I>-derivable approach to an effec- 
tive three-dimensional field theory for the zero-frequency 
mode that has been obtained by dimensional reduc- 
tion. The results are comparable to those obtained in 
SPT [13 and the ^-derivable approximation in 3+1 di- 
mensions |3J|. 

In this paper, we solve the three-loop ^-derivable 
approximation for QED. This is done by applying the 
strategy developed in Ref. |3J| . It consists of expanding 
the sum-integrals systematically in powers of e and m/T 
where ?77, is a variational mass parameter of order eT. 
Our result for the pressure reproduces the weak-coupling 
expansion to order g^ and is gauge invariant. 

The paper is organized as follows. In Sec. II, we 
briefly discuss the application of the <I>-derivable ap- 
proach to QED and the general framework developed 
in Ref. ^^ to solve it systematically. In Sec. HI, we 
solve the two-loop 4>-derivable approximation and dis- 
cuss the issue of gauge dependence. In Sec. IV, we solve 
the three-loop ^-derivable approximation. We summa- 
rize and draw some conclusions in Sec. V. There are two 
appendices where our notation and conventions are given 
and where we list the sum-integrals and integrals that are 
needed. 



II. ^-DERIVABLE APPROXIMATIONS 

In this section, we briefly discuss the 2PI effective 
action formalism and ^-derivable approximations. 
The Euclidean Lagrangian of massless QED is 



1 



^ = T^M. + V'Tm-D/.V' + -Cgf 



(1) 



where F^u — d^A^ — d^A^j^ is the field strength tensor, 
D^i = dfj^ + ieA^ is the covariant derivative, and e is 
the electric coupling. £gf is the gauge-fixing part of the 
Lagrangian. In general covariant gauge, the gauge-fixing 



part of the Lagrangian is 



^gf = 21^^'^^'^)^ 



(2) 



In the remainder of this Sec. and in Sec. Ill, we keep 
^ general in order to discuss the problem of gauge de- 
pendence. In Sec. IV, we specialize to Feynman gauge 
(^ = 1), which by far is the easiest gauge for practical 
calculations. 

The thermodynamic potential fi of QED is 

n[A,S] = ixrlogA-i-TrlogS-i-TrlogA-i^i 



-TrHA^ 
2 

-*[A,5], 



TrSS' + Trn.hAgh 



(3) 



where A^^(P) and S{P) is the exact photon and electron 
propagator, respectively, and Agh(P) is the propagator 
for the ghost. np^(P) is the polarization tensor and S(P) 
is the electron self-energy. We can then write 

A;J(f) = [A%{P)r' + u,^iP), (4) 

S-\P) =f + S(P), (5) 

where A'^ (P) is the free propagator in covariant gauge: 
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FIG. 3: One- and two-loop photon self-energy graphs. 
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FIG. 4: One- and two-loop electron self-energy graphs. 



The exact propagators satisfy the variational equa- 



tions 
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(7) 
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Using Eq. Q, the variational equations Q and (|SJ) can 
be written as 



n^.(P) = 2 



s(P) = - 



(5$[A,S'] 
5A,,,{P) ' 
J$[A,6'] 

ss{p) ■ 



(9) 
(10) 



A° (P) 



P2" 



(1 
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The trace in Eq. Q is over Dirac and Lorentz indices 
as well as space-time. In covariant gauges, the ghost 
field decouples from the other fields and so the ghost 
self-energy ngh(P) vanishes identically ^ . The functional 
<I>[A, S] is the sum of all two-particle irreducible vacuum 
diagrams. We define the n-loop ^-derivable approxima- 
tion f2„ to the thermodynamic potential O as the trun- 
cation of the action functional after n loops. The two- 
particle irreducible vacuum diagrams are shown diagram- 
matically in Fig. |21 up to three-loop order. The corre- 
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FIG. 2: ^-derivable two- and three-loop skeleton graphs. 

sponding self-energies that are obtained by cutting a line, 
are shown in Figs. 13 and 01 




^ In nonabelian gauge theories, the ghost does not decouple in 
covariant gauges which makes the calculation significantly more 
involved. In Ref. |32| the authors are employing the temporal 
axial gauge in which the ghost does decouple. Howeve r, ther e 
are other problems with this gauge at finite temperature l39tl4Cl . 



In QED, we know that thermal fluctuations generate 
a mass m for the zeroth component of the gauge field Aq 
which is of order eT and screens the interactions. The 
strategy for solving the n-loop ^-derivable approximation 
is to introduce a mass variable which is of order eT and 
then calculate the sum-integrals as double expansions in 
e^ and m/T. This strategy was developed in Ref. |3J| in 
order to solve the three-loop ^-derivable approximation 
in scalar field theory. It turns out that the gap equa- 
tions and (fTU)) have a recursive structure that allows 
us to solve for their dependence of the external momen- 
tum P. We follow Braaten and Petitgirard and choose 
the Debye mass as the mass parameter. The Debye mass 
is the solution to the equation 



P 



noo(o,p) 



0, 



p 



(11) 



In the variational equations ^ and H10() . there are two 
important mass scales. One is soft and is of order eT. 
This scale is set by the Debye mass m. The other is the 
hard scale of order 27rT and is set by the nonzero Mat- 
subara modes. We will assume the coupling is sufficently 
small so that the scales m and 27rT are well separated. 
This allows one to expand the sum-integrals in powers of 
e^ and m/T . The gap equations will then be solved in the 
two momentum regions separately. For hard momentum 
P, we expand the polarization tensor as follows: 



n^.(P) = 



+nf,'J(P)- 



°-* re(P) 



(12) 



where nj^j7(P) is of order T^{m/Ty. Similarly, the elec- 
tron propagator is expanded as 



, (13) 

where E"''^(P) is of order T'^{ra/T)^ . For soft momentum 
P ~ (0, p), we expand the longitudinal part of polariza- 
tion tensor as follows ^: 



noo(0,p) = m^ + e^ 



\v) 



-e^ ^4.0^ 



\P) + -] 
where cr'^'^{p) is of order m?{mlT)^ . 

For hard momentum, we can expand ^Trlog A^^ 
about the free propagator, since the self-energy is per- 
turbative corrections starting at order e^. This yields 



^TrlogA-i = lid+l)^logP' + ^e^^ 



p2 



H 



p2 



in2.o(p)n2^o(p) 



pi 



(15) 



The gauge-dependent terms in Eq. (|15|l drop out since 
the photon self-energy in QED is transverse to all orders: 

(16) 
,(17) 



p^n^AP) = 0. 

For soft momentum, the expansion is 



-TrlogA^i 
2 ^ 



¥ 



2 



log (p2 



log [p 
1 



-e^T 



+ 6^0-2, o(p) 

'^2fi{p) , 



,p - jpP^ + m^ 

Again we have used the transversality of the photon prop- 
agator to eliminate the gauge-dependent terms. We do 
not need the expansion of TrllZ? since we will use the 
gap equation to eliminate this term. 

We also need the expansions for Trlog5~^ and 
Tr ES*. Since the electron momentum is always hard, 
we can expand about the free propagator and obtain 



Trlog5"^ = 2 



logp2 



{P} 

2 



4 

J{P} 



rv2.o 



Tr 



I]2^0(P)f E24(P)|& 



i 

T{p} 



p2 

S2."(P)fI]2'0(P)f 



p2 



TrSS" = 



Tr 



P} 



pi 

E2-0 



S2'l(P)f 



p2 



p2 



^ r s^-o(p)fs^.o(p)f 



j{p} 



pi 



(18) 



, (19) 



^ Since the infrared limit of the other components of 11^^(0, p) 
vanishes, the corresponding contribution to the free energy also 
vanishes. 



where the trace on the right-hand side is only over Dirac 
indices. 

The contribution from the ghost field is as usual 

(20) 



TrlogAgi^i-^logp2 



By inserting the expansions for the self-energies into 
the gap equations and expanding systematically in pow- 
ers of e and m/T, we obtain expressions for nj]^'^(P), 
(14') S"''^(P) and <T^'i^{p). By matching coefficients of e" on 
both sides and solving the equations simultaneously and 
recursively. 



III. TWO LOOPS 



In the two-loop <i>-derivable approximation, there is 
only a single diagram contributing to ^[D,S] which is 
the left diagram in Fig. |21 The two-loop thermodynamic 
potential fl2 is 



n2[A,S] = -TrlogA^i -TrlogS'^^-TrlogA 



-1 

gh 



-TrBA + TrSS'+ie^ 



P{Q} 



xTr[5(Q)7^5(P + Q)7'^A^,(P)] . (21) 

The gap equations are obtained by varying the thermo- 
dynamic potential il2 with respect to np^(P) and S(P): 



nM.(i^) 



i 

J{Q} 



TT[S{Q)rSiP + Q)r] , (22) 



S(P) = e'jf j'^S{P + Qh''A,,,{Q) . (23) 

It follows from the coupled gap equations that both 
n^!v(P) and E(P) are nontrivial functions of the external 
momentum P. 

The gap equation (|22|l can be used to simplify equa- 
tion (EU for ^2' 



^2 = iTrlogA-i-TrlogAgh'-Trlog5-i 
-t-TrSS". 



(24) 



Substituting the expansions for the various terms 
into 1)24(1 and truncating at the appropriate order, we 
obtain 

^2 = ^{d-l)^\ogP' + ^Tj\og{p' + m') 

1 _r u'^'"(P) 



-2jl logp2 + ie2^ 
J{P} ^ JP 



p2 



(25) 



We note that the function Y?'^(P) drops out. In the two- 
loop ^-derivable approximation, we only need the trace 



of n^'^(P), while in the three-loop ^-derivable approxi- 
mation, we need the function itself. 

The solution to the gap equations for hard momen- 
tum to order e^ are obtained by using bare propagators 
in the loops. In this manner, we find 



^t'iP) 



{Q} 
til 



8QaQi, 



2P^d., 



_ 4J^ 
, 4P^Q. + 4P,Q^ 



QHP + Q)' 



Q^P + Q? 



(26) 



We also need to solve the gap equation for soft mo- 
mentum in order to determine the Debye mass. Through 
order e'^, the longitudinal part of polarization tensor at 
zero frequency reads 



noo(o,p) 

where the Debye mass is 
167r2 



m 



l)e' 



{Q} 



aT^ 



(28) 



where a = e^/(47r)^. The thermodynamic potential 
through order e^ then reduces to 



^2 = ^(d-1 



\ogP^ + -T j log(p2 

2^ \ogP^-[d-l)e^ 
J{P\ 

JF 



IP{Q} 



P^Q^ 



{PQ} 



P^Q^ 



(29) 



Using the expressions for the integrals and sum-integrals 
in the appendices, this reduces to 



2rp4 



n. 



180 



50 
— ( 

11 



320%/3 
33 



-a 



3/2 



(30) 



Eq. (|3U|I agrees with the weak-coupling result through 
order e^ (6j. Thus the two- loop <i>-derivable approxima- 
tion sums up the leading contribution from the plasmon 
diagrams. 

We close this section by discussing the problem of 
gauge dependence that arises when going beyond order 
e^. For example, to calculate the Debye mass to order e'', 
we need to include the function E^'°(P) in the dressed 
electron propagator on the right-hand side of the gap 
equation H22|) . This function is 



E2,0(p) = (l-d)^ 



-m 



gQ^P + Qf 



-(1-0 



Q\P + Qf 



.(31) 



The function E^^''(P) arises from hard photon momenta 
in the self-energy graph in Fig.0] We note that it is gauge 
dependent and this is due the photon line in the one-loop 
self-energy graph shown in Fig.^ Since E^'°(P) is gauge 
dependent, this introduces a gauge dependence at order 
e^ in the longitudinal part of polarization tensor. At zero 
frequency, one finds: 



noo(o,p) 



-4((i 



1)^4 



T{Q} Q' 



+ -{d-l)e 



J{Q} y 



+ |(d-l)(d-3)e4^ 



T{Q} Q' 



% i?2 9|il> R' 



(27) +8(l-e)e^rf [ 



p2 



{R} 

2Pi 



{p}Q iP'^Q^P + Q? P'Q'l ■ 

(32) 

The equation for noo(0, p) is ultraviolet divergent and re- 
quires renormalization. The divergence proportional to 
p^ is removed by wave-function renormalization in the 
usual manner (see also Sec. IIVJI . The other divergence 
which arises from the last line in Eq. 1)32(1 can only be re- 
moved by a gauge-dependent renormalization of the cou- 
pling constant '^. This would lead to a gauge-dependent 
gap equation and Debye mass. Eq. ((32|l depends on the 
gauge-fixing parameter, but this is not in contradiction 
with the fact that the photon self-energy is manifestly 
gauge invariant. The point is that in the two-loop <&- 
derivable approximation, we are not including all contri- 
butions to n^i^(P) of order e^. This is done when one 
considers the three-loop approximation and we have ex- 
plicitly checked that the gauge dependence cancels alge- 
braically as we include the two-loop self-energy graph in 
Fig. El 



IV. THREE LOOPS 

The three-loop ^-derivable approximation to the 
free energy is 

n,[A,S] = iTrlogA-i-Trlog5-i-TrlogA;i^i 



gh 



1 

2 



TrHA + TrES* 

e-2l TT[S{Q)rSiP + Qh''A^,iP)] 
^ JP{Q} 

]e^jl Tr [S{Q)rS{R)rS{R - P) 

4 Jp(OR.\ 



^ Only in Feynman gauge is the polarization tensor finite after 
wave-function renormalization. 



xj^SiQ ~ Ph^A^,{P)A^p{Q - R)] . 

(33) 

The gap equations are again obtained by varying H33|l 
with respect to the photon and electron self-energies: 






'{Q} 



xS{R - Ph^SiQ - P)7^] Ao,i3{Q - R) , 

(34) 

S(P) = e'j{ rS{P + Q)rA^,{Q) 

JQ 

rS{R)rS{R-Q)i'' 
Q{R} 

xS{P-Q)-i^A^p{Q)A^,{R-P). (35) 

The gap equation for 11^^ (P) can be used to simplify the 
expression for the thermodynamic potential: 

n^[A,S] = iTrlogA-i-Trlog^-i-TrlogAgi^i 
+TrSS' 
-ie^rf Tr [S{Q)YS{R)rS{R - P) 

4 JP{QR} 

xYS{Q - P)^^A^,{P)A^p{Q - R)] .(36) 

Substituting the expansions p5|l -H19 |l into l|36(l . we ob- 
tain 



!^3 = 



i(d - 1)^ logp2 + iry log(p2 + jn") 



2 JpP^ + rn? 



2jf logp2 
T{P} 



H 



^m 



2" fp P2 



n 



n^'^°(P) , n4-,i(p) 



P2 



-K 



P2 



1^ n^°(p)n^-°(p) 



4^rf Tr 
2 T{P} 



P4 

I]2'°(P)fS2'0(P)f 

P4 



-H2 



S2,0(p)|Ss2.1(P)f I]2.1(P)fE24(p)|& 



p4 



p4 
1 



4(rf-l)(5-d)e^^^^^^^^^,^^^Q^^^, 



+ {d~l){d^3)e^j{ 



+8{d-l) 






%q)rQ'R''{Q + R)^ 



pP^ + m? 



(37) 



In the three-loop ^-derivable approximation, we only 
need the trace of the functions n^'°(P) and n^'J(P). 
These traces are significantly simpler to calculate than 
the functions themselves. Furthermore, it turns out that 
the order-e^ term in equation (|35|1 for the electron self- 
energy first contributes to the pressure at order e^ and 
so it is not needed. In addition to the functions listed in 
Sec. Ill, we need the following functions for hard P 



I]2.i(P) = 



2po7o- f 

P2 



T 



1 



g2 + rn? 



(38) 



J{Q}R 



1 



p2 



Q4i?2 Q4p2(Q + p)2 

1 2P-P 



-4(d-l 



{QR} 



P2 



Q4i?2(p + Q)2 Q4p2 

4(Q ■ i?)2 

Q2p2(^_p)2(Q_p)2(g_^)2 

P4 

Q2p2(p_p)2(Q_p)2(Q_^)2 

4P2 



{Qfl} 



Q2p2(Q„p)2(Q_p)2 

d-7 P2 

2~Q2p2(p_p)2(Q_p)2 

^^'^ " ^^P2(g_p)2(Q_p)2 

2Q2 

i?2(P-P)2(Q_p)2(Q_p)2 

2 



i?2(P-P)2(Q_p)2 

" 2P2 



(39) 



n^'^i(P) = 2(d-l)5f 



8P^go' 



QHP + Q? Q^ 

2„2 



^pmoP'^ 



M 

Q6 Q6(p + g)2 Q4(p + Q)4 



4p2g2 



Q4(p + g)'^ 



T 



1 



r^ + ?7i2 



(40) 



The functions S2^i(P) and n^^J(P) arise when the pho- 
ton momentum in the relevant Feynman diagram is soft, 
while nl^'^(P) is when the photon momentum is hard. 

Through order e^ , the longitudinal part of polariza- 
tion tensor at zero frequency can be written as 



noo(o,p) 






+Aid-l){d-3)e'j{ ij 
T{Q} Q 



% K' %} R' 



+4(d-l)(d-3)e 



'±J. 



T{P} 



(41) 



We have explicitly checked that (|41|l is independent 
of the gauge-fixing parameter i^ in contrast to the gauge- 
dependent expression H32II . Eq 1411) can be rewritten as 



noo(0,p) - m2 + eV'°(p), 
where the Debye mass satisfies 

1 



(42) 



-4(d- l)e' 



{Q} 



Q^ 






+Aid-l)id-3)e'jf ij 
T{Q} y 



+4(.-l)(.-W i,(T/-^ 



and 



'^^''^(p) 



^(d-l)e2(p2+TO' 



T{Q} y 



(44) 



The gap equation (|43|l is ultraviolet divergent and re- 
quires renormalization. The divergence is proportional 
to p^ and is removed by wave-function renormalization: 



Za = l- 



Uir^e 



(45) 



After having renormalized the static polarization tensor 
noo(0,p), the gap equation ((TT|) reduces to 



167r2 „ 8 

T^a< 1- - 

3 I 3 



log^+7 + 21og2 + - 



-18 



f— ) 



(46) 



The result for the Debye mass (|46|) agrees with the weak- 
coupling result 7, 41] through order e^. 

Inserting the expressions (|38|I - H4U() and 144|) 
into H37|l . the three-loop ^-derivable approximation then 
becomes 

«3 - ^id-l)^\ogP' + ^Tj\og{p' + m') 

-2jl logp2_e2(d-l) 
T{p} 



-^'-'Hn^^ 



T|n> Q2 % Q2 



'{Q} 
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-'''-'H,,.,P^Q^R' 



'2{d~l)e 

--{d^~8c 
2^ 



■t 



1 



11)6^ 

1 
f{PQR} P'Q^BHP + Q + W 

-2(d-i)Vi ^^ 



ny 



4 



-16e' 



'il 



iQR? 



rp{QB.} P^Q^RHP + QYiP + R? 



^Yf 1- 



-{d-l){d-3)e-^^ —JT 



1 



q^ + m^ 



(47) 



(43) Note that we have kept a term that is proportional to 
e^m? and first contributes at order e^. This term arises 
from three- loop diagrams where both photons are soft ■*. 
This contribution is manifestly gauge invariant and we 
include this selective resummation in our final result. It 
is interesting to note that the only contribution at order 
e^ comes from the Debye mass; all the other contributions 
cancel algebraically 

The expression for fia is ultraviolet divergent. The 
divergences can be eliminated by renormalizing the cou- 
pling constant. This is done by the substitution e^ — > 
ZgE^, where 



1 



127r2e 



(48) 



Using the expressions for the sum-integrals and integrals 
listed in the appendices, we obtain 
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(49) 



Using the expression for the Debye mass in Eq. (|46|l . one 
can show that the three-loop ^-derivable approximation 



* The complete e'^m? contribution can also be obtained in a two- 
loop calculation using the dimensionally reduced theory of QED 



(electrostatic QED) derived in f^lill 



agrees with the weak-couphng expansion through order 

The renormahzation group equation that follows 
from (gHJ is 



M' 






e 
6^ 



(50) 



which coincides with the standard one-loop running of 
the coupling. 

In Fig. 13 we show the two and three-loop $- 
derivable approximations to the pressure normalized to 
that of an ideal gas shown as dashed and solid lines, re- 
spectively. In the three-loop approximation, the band is 
obtained in the usual manner by varying the renormal- 
ization scale fi. The three-loop band is slightly narrower 
when compared to the e^-band in Fig. ^ The approx- 
imations also seem to be slightly more stable than the 
successive weak-coupling approximations. However, the 
final result does not seem to be a dramatic improvement 
over the e^ result. 




lancu and Rebhan in the case of QCD |30|, |3l|, |33| ■ In 

scalar field theory, the coupli ng is running incorrectly by 
a factor of three [Sfl 113, IHj Ull • In the three- loop <I>- 
derivable approximation, the expressions for the Debye 
mass and the free energy require wave-function renor- 
mahzation and renormalization of the coupling constant, 
respectively. The running of the resulting renormalized 
coupling agrees with the standard one-loop running in 
QED. 

We have also considered the problem of gauge depen- 
dence within the 2PI effective action formalism. We gave 
an explicit example of the how gauge dependence arises 
in the one-loop gap equation for the photon propagator 
when one truncates the gap equations at order e'*. Our 
calculation is in agreement with general results on the 
gauge dependence of ^-derivable approximations [33 ■ 

The method could also be used to solve the three- 
loop ^-derivable approximation in QCD with an accu- 
racy of order g^; however, given that the final three- 
loop ^-derivable result in QED does not seem to dra- 
matically improve the scale variation at large coupling it 
is questionable whether this would be a worthwhile en- 
deavor. Going beyond three- loops in scalar theory as well 
as gauge theories would be very difficult. One problem 
is that there are new four-loop sum-integrals of order g^ 
that have not yet been evaluated. In nonabelian gauge 
theories, there is the additional problem that the free 
energy at this order is sensitive to the nonperturbative 
momentum scale g^T which is associated with screening 
of static magnetic fields. This may cause the expansion of 
the three-loop <I>-derivable approximation to break down 
beyond order g^. 
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FIG. 5: Two- and three- loop ^-derivable approximations to 
the pressure normalized to that of an ideal gas as a function 
of e(27rT) shown as dashed and solid lines, respectively. The 
three-loop band is obtained by varying the renormalization 
scale, /i, by a factor of two around the central value fi — 2nT. 
Note that the scale is different than in Fig. 
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APPENDIX A: SUM-INTEGRALS 



V. SUMMARY 

In this paper we have solved the three-loop $- 
derivable approximation for the Debye mass and the free 
energy in QED by systematically expanding the sum- 
integrals in powers of e^ and m/T. The results are 
accurate to order e^. In the two- loop ^-derivable ap- 
proximation, both the thermodynamic potential and gap 
equation are finite, and so there is no running of the cou- 
pling. The solution to the gap equation is trivial; the De- 
bye mass is given by its weak-coupling expression. This 
was also the case in the two-loop calculation of Blaizot, 



In the imaginary-time formalism for thermal field 
theory, the 4- momentum P — (PoiP)is Euclidean with 



p2 



p2 



p . The Euclidean energy po has discrete 



values: Pq = 2n7rT for bosons and Pq = {2n + 1)'kT for 
fermions, where n is an integer. Loop diagrams involve 
sums over Pq and integrals over p. With dimensional 
regularization, the integral is generalized to d = 3 — 2e 
spatial dimensions. We define the dimensionally regular- 
ized sum-integral by 
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where 3 — 2e is the dimension of space and /i is an 
arbitrary momentum scale. The factor {e"' /AttY is in- 
troduced so that, after minimal subtraction of the poles 
in e due to ultraviolet divergences, fi coincides with the 
renormalization scale of the MS renormalization scheme. 



1. One- loop sum-integrals 

The specific one-loop sum-integrals needed are 
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The errors are all one order higher in e than the 
smallest term shown. The calculations of these sum- 
integral is standard. 



3. Three-loop sum-integrals 

The three-loop diagrams needed are 
1 
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2. Two- loop sum-integrals 



The two-loop sum-integrals that are needed all van- 



ish: 
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pqPW{P + W 
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JF 
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The errors are all of order e. Details of the calcula- 
tion of these two-loop sum-integrals can be found in e.g. 
Ref. 0. 



The errors are all of order e. The calculation of these 
three-loop sum-integrals was done in Ref. pj and details 
can be found there. 



APPENDIX B: INTEGRALS 

Dimensional regularization can be used to regularize 
both the ultraviolet divergences and infrared divergences 
in 3-dimensional integrals over momenta. The spatial di- 
mension is generalized to d = 3 — 2e dimensions. Integrals 
are evaluated at a value of d for which they converge and 
then analytically continued to d = 3. We use the inte- 



gration measure 



An 



(27r)3 



(Bl) 



We require a few integrals in that appear in the soft 
sector. The momentum scale in these integrals is set by 
the Debye mass m. The one-loop integrals needed are: 



log (p2 



m 



(B2) 



m 

477 



(B3) 



The errors are all of order e. The calculation of these 
integrals is standard. 
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